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FOREWORD 


The  work  described  herein  was  authorized  under  Task 
4C04-15-029 “09,  Mathematical  Research  on  Dissemination  (U).  The  work 
was  started  in  September  1958  and  completed  September  I96I. 

Notice 

Reproduction  of  this  document  in  whole  or  part  is  prohibited 
except  with  permission  of  the  issuing  office;  however,  ASTIA  is  authorized 
to  reproduce  the  document  for  U.  S.  Governmental  purposes. 

Disposition 

When  this  document  has  served  its  purpose,  DESTROY  it.  DO 
NOT  return  the  document  to  U.  S.  Army  Chemical  Research  and 
Development  Laboratories. 


DIGEST 


Slide  rules  have  been  contrived  for  evaluating  O.  G.  Sutton's 
mathematical  models  for  estimating  the  downwind  travel  and  area  coverage 
from  an  instantaneous  point  source  of  vapor  or  very  fine  aerosol 
disseminated  and  sampled  at  ground  level. 

Using  these  slide  rules,  the  nonmathematician  can  employ  the 
mathematical  theory  for  making  preliminary  weapons -effects  predictions 
necessary  in  handling  everyday  problems  encountered  in  research  and 
development  of  chemical -weapons  systems. 
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CLOUD  TRAVEL  SLIDE  RULES  FOR  THE  DETERMINATION 
OF  DOWNWIND  TRAVEL  AND  AREA  COVERAGE 


I.  INTRODUCTION. 

In  the  research  and  development  of  a  chemical -weapons  systems, 
determination  of  the  weapons  characteristics  for  producing  the  optimum 
weapons  effects  is  basic.  The  large  number  of  meteorological  and  physical 
conditions  to  be  considered  in  a  weapons -evaluation  study  precludes  field 
experimentation  as  an  adequate  technique  for  evaluation.  Thus,  mathematical 
models  have  been  developed  to  estimate  expected  results;  however,  the 
mathematical  forms  of  these  models  are  such  that  considerable  numerical 
manipulations  are  required  in  determining  a  set  of  data.  Thus,  the  mathe¬ 
matical  models  for  estimating  downwind  cloud  travel  and  area  coverage  have 
been  reduced  to  the  form  of  a  mechanical  slide  rule  in  order  to  simplify 
sufficiently  some  aspects  of  weapons -effects  predictions  in  such  a  manner  that 
the  nonmathematician  can  apply  the  theory  to  everyday  problems  encountered 
in  research  and  development  of  chemical -weapons  systems. 

II.  THEORY. 

A.  Diffusion  Equations . 

Sutton's*  nonisotropic  diffusion  equation  for  predicting  total  dosage 
downwind  from  an  instantaneous  point  source  composed  of  vapor  or  very  fine 
aerosols  disseminated  at  ground  level  is 


where 


D(x,y,z) 


2Q 


-  _  _  2-n 
TT  CyC^UX 


exp 


(1) 


D  =  dosage,  mg  min/cu  m 
Q  =  total  mass  of  agent  dispersed,  mg 


O.  G.  Sutton,  The  Problem  of  Diffusion  in  the  Lower  Atmosphere. 
Quart.  J.  Roy.  Meteorol.  Soc.  73,  268  (1947), 


i 
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Preceding  Page  Blank 


generalized  coefficients  of  diffusion  crosswind  and 
vertically,  m 

mean  wind  speed,  m/min 

distance  downwind  from  source,  m 

distance  crosswind  from  mean  wind  track  or  x-axis , 

sam  'ng  height  above  ground,  m 

parameter  related  to  atmospheric  stability 


It  is  seen  from  equation  (1)  that  the  ground-level  dosage  is 

given  by 


D(x,y,o)  = 


2Q 


TrC^C„ux^-'^ 

y  ^ 


exp 


-y 


C^x^-n 


The  expression  describing  the  contour  of  dosage  D  is  found 
by  rearranging  equation  (2) 


y 


CyX 


(2-n) 

~Y 


1 

2 


where 


,  Sutton's  expression  for 


i{2-n) 

Since  y  =  o  at  x  =  o  and  at  x  =  a^ 
the  area  enclosed  by  this  contour  is  given  by 


1  r-r  r  2Q 

(2.n)  ^  y  [TTCyC2nD(x 

A  =  2  r  ydx  =  - 


2Q  1 
TTCyC2nD(x,  y.  o)J 


(4-n) 

(4-2n) 


From  equation  (2),  it  is  seen  that  the  ground-level  dosage  along  the 
mean  wind  track  is 


D(x,  o, o) 


itC„C„ux^-^ 


The  ground -level  semiwidth,  y,  of  the  cloud  at  axial  distance  x 
downwind  from  the  source  is  the  distance  from  the  x-axis  at  which  the  dosage, 
D(x,y), is  the  fraction  K  of  the  axial  value,  D(x,o).  With  equations  (2)  and  (5), 
it  is  easily  deduced  that 

(2  -n) 

y  =  CyX  ^  \y”-lnl?  (6) 


!3.  Diffusion  Coefficients. 


Sutton's*^  definitive  expressions  for  the  diffusion  coefficients  Cz 
and  Cy  for  airflow  over  an  aerodynamic  ally  smooth  surface  are 


r<2  . 


(l.n)(2.n)a"  \u^ 


\  1-n 


'W 


*  O.  G.  Sutton.  The  Theoretical  Distribution  of  Airborne  Pollution  From 
Factory  Chimneys.  Quart.  J.  Roy,  Meteorol.  Soc.  73,  426-436  (1947). 


where 


V  =  kinematic  viscosity  of  air  (equal  to  0.  143  x  10*^  sq  m/sec  at 
1,000  mb  and  10 “C)  ... 

u  =  mean  wind  speed,  m/sec 

w*,  V*  =  eddy  velocities  in  vertical  and  crosswind,  respectively. 

n  =  a  dimensionless  stability  parameter. 

Sutton  assigned  values  of  n,  U,  and  the  ratios  (^ ' /IT)  i 

as  being  representative  for  each  of  four  stability  conditions  shown  in 
table  1. 

TABLE  1 

VALUES  OF  n,  u,  AND  w '/ll  FOR  VARIOUS  STABILITY  CONDITIONS 

(Sutton) 


Stability  condition 

n 

H 

(w'/u)n/(3V’/u)i. 

4  . 

m/sec 

Large  lapse 

1/5 

7 

2 

Neutral 

^  1/4 

5 

1 

Moderate  inversion 

1/3 

3 

1/2 

Large  inversion 

1/2 

2 

1/4 

For  practical  purposes,  the  following  (table  2)  may  be  used  as 
a  guide  for  defining  stability  conditions. 

Sutton  has  also  assigned  values  for  (table  3)  at  various  heights 
ac  being  appropriate  for  neutral  conditions  in  a  wind  of  5  m/sec  over  smooth 
terrain. 
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TABLE  2 


STABILITY  PARAMETERS 


Temperature  gradient 
T4.  0  m  -  Tq.  5  nri 

Stability  conditions 

n 

•F 

AT  <-0.4 

Lapse 

0.20 

-0.4  <  AT  <0.4 

Neutral 

0.25 

0.4  <  AT  <2.0 

Moderate  inversion 

0.  33 

AT  >2.  0 

Large  inversion 

0.50 

TABLE  3 

VALUES  OF  C,  FOR  VARIOUS  HEIGHTS  WHEN 
- SL_ - ti:; - 

n  =  =  5  m/sec 

(Sutton) 


Height 

m 

n 

0-25 

0.  12 

50 

0.  10 

75  1 

0.09 

100 

0.  07 
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Sutton*  assumes  Cy  =  at  a  height  of  E5  m  and  above.  In  order 
to  obtain  values  of  Cy  for  the  various  stability  conditions  at  heights  <25  m, 

however;  values  of  the  ratios  (v '/Tr")ji/(v  •/'u)^  are  needed.  Best**  con¬ 
cluded,  from  data  collected  from  extensive  tests,  that  the  ratio  of  lateral  to 
vertical  gustiness,  v'/w’,  is  1.81  irrespective  of  wind  velocity  and  temper¬ 
ature  gradient.  Using  this  result  with  equation  (7),  the  following  relationship 
is  obtained: 


Cy  = 


(8) 


With  equations  (7)  and  (8)  and  tables  1  and  3,  values  of  the 
diffusion  coefficients  may  be  determined  for  the  appropriate  conditions  as 
listed  and  specified  by  Sutton. 


By  assuming  that  the  ratios  (w'/u)  and  (v'/H)  are  constant  for 
various  values  of  u  (Suttonf  says  they  are  "practically  independent  of  the 
magnitude  of  H. "),  extensions  of  these  coefficients  to  include  other  wind 
speeds  may  easily  be  accomplished.  For  constant  n,  it  may  be  deduced 
from  equation  (21)  that 


(9) 


The  above  outlined  procedure  was  used  to  determine  the  diffusion 
coefficients  used  here.  A  more  detailed  explanation  of  the  procedure  along 
with  tabulated  parameter  values  for  a  variety  of  conditions  are  contained  in  a 
previously  published  report,! 


*  O,  G.  Sutton,  The  Theoretical  Distribution  of  Airborne  Pollution  From 
Factory  Chimneys.  Quart.  J.  Roy.  Meteorol.  Soc.  73^,  426-436  (1947). 

**  A.  C.  Best.  Transfer  of  Heat  and  Momentum  in  the  Lowest  Layers  of 
the  Atmosphere.  Geophys,  Memoirs  No.  65.  Air  Ministry  Meteorol, 
Office,  London.  1935. 

t  O.  G,  Sutton.  The  Problem  of  Diffasion  in  the  Lower  Atmosphere. 
Quart.  J.  Roy,  Meteorol.  Soc.  73,  268  (1947). 

t  D,  ,0.  ..Egner  ^  CRDLR  3002.  A  Preliminary  Theoretical  Study 
■  of  tne  Cloud  Travel  of  Aerosol  Particles  Having  Diameters  Between 
20  and  180  Microns.  July  I960, 
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III.  DESCRIPTION  AND  USE  OF  CLOUD  TRAVEL  SLIDE  RULES. 


A.  Dosage  Field  Slide  Rule. 

The  dosage  field  slide  rule  is  used  to  determine  total  dosage  at 
any  point,  x,  y,  downwind  from  an  instantaneous  point  source  disseminated 
at  ground  level  over  smooth  terrain  under  various  wind  speeds  and  atmos¬ 
pheric  stability  conditions.  The  front  side  of  the  slide  rule  is  used  to 
determine  the  total  dosage  along  the  mean  wind  track,  D(x,  o),  and  the  reverse 
side  is  used  to  determine  the  semiwidth,  y,  of  the  cloud,  containing  a  desired 
dosage  at  a  specified  downwind  distance.  Thus,  the  entire  dosage  field  may 
be  established  by  determining  the  dosage  values  at  a  sufficient  number  of 
points.  The  scales*  for  the  dosage  field  slide  rule  are  shown  in  figure  1. 


The  following  two  examples  illustrate  the  use  of  both  sides  of 
this  slide  rule.  Given  any  four  parameters  included  on  a  side,  the  fifth 
parameter  value  may  be  obtained  by  rearrangement  of  the  steps  outlined 
in  the  examples. 

PROBLEM  1:  Find  the  ground  level  total  dosage  along  the  mean  wind 

track  at  lOOm  downwind  from  al,000  gm  instantaneous 
point  source  disseminated  at  ground  level  under  the  follow¬ 
ing  conditions:  Wind  speed  at  2  m  height  =  2  m/sec, 
neutral  stability  (n  =  0,25). 


SOLUTION:  1.  Use  front  side  of  dosage  field  slide  rule. 

2.  Set  hairline  of  cursor  on  Q  value  of  1000  on  D  scale. 

3,  Move  center  slide  until  assumed  stability  parameter,  n,  of 
0.25  on  cursor  intercepts  assumed  wind  speed,  H,  of  2  m/sec 
on  C  scale. 


4.  Move  cursor  until  n  value  of  0.25  on  cursor  intercepts 
desired  downwind  distance,  x,  of  lOO  m  on  B  scale. 


5.  Read  dosage  on  A  scale  indicated  by  hairline  of  cursor 
(D  =  65  mg  min/cu  m). 

PROBLEM  2;  For  the  same  source  strength,  wind  speed,  and  stability 
condition  assumed  in  problem  1,  find  the  semiwidth  of  the 
cloud,  y,  containing  a  dosage  of  6.5  mg  min/cu  m,  at  a 
downwind  distance  of  100  m. 


*  The  unlabeled  values  on  the  cursors  are  values  of  the  stability  parameter, 
n,  as  given  in  table  2. 
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DOSAGE  FIELD  SLIDE  RULE 


SOLUTION:  1.  Use  back  of  dosage  field  slide  rule. 

2.  In  problenn  1,  the  center-line  dosage  at  100  m  downwind 
was  found  to  be  65  mg  min/cu  m.  Thus, 


K  = 


D(x,  y)  _  6^ 
D(x,  o)  “  65 


=  0.1 


3.  Set  hairline  of  cursor  at  K  =  0.  1  on  E  scale. 


4.  Move  center  slide  until  assumed  stability  parameter,  n, 
of  0.25  on  cursor  intercepts  assumed  wind  speed,  11,  of 
2  m/sec  on  C  scale. 

5.  Move  cursor  until  n  value  of  0.25  on  cursor  intercepts 
desired  downwind  distance,  x,  of  100  m  on  B  scale. 

6.  Read  semiwidth  of  cloud  on  F  scale  indicated  by  hairline 
of  cursor  (y  =  17  m). 


B,  Area  Dosage  Slide  Rule. 

The  front  side  of  the  area  dosage  slide  rule  is  used  to  determine 
expected  area  coverage  to  a  desired  dosage  from  an  instantaneous  po^nt 
source  under  the  same  range  of  wind  speeds  and  stability  conditions  con¬ 
sidered  by  the  dosage  field  slide  rule.  The  reverse  side  is  an  exact  duplicate 
of  the  front  side  of  the  dosage  field  slide  rule  and  requires  no  further 
explanation.  The  scales*  for  the  area  dosage  slide  rule  are  shown  in 
figure  2. 


The  following  example  illustrates  the  use  of  this  slide  rule  for 
determining  area  coverage.  As  in  the  instance  of  the  dosage  field  slide 
rule,  given  any  four  of  the  parameters  included,  the  fifth  parameter  value 
may  be  obtained  by  rearrangement  of  the  steps  outlined  in  this  example. 

PROBLEM:  Find  the  area  covered  to  a  dosage  of  lOO  from 

a  1,  000  gm  instantaneous  point  source  disseminated  at 
ground  level  under  the  following  conditions: 

Wind  speed  at  2  m  height  =  2  m/sec 

Neutral  stability  (n  =  0,25) 


*  The  unlabeled  values  on  the  cursors  are  values  of  the  stability  parameter, 
n,  as  given  in  table  2. 
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FRONT 


SOLUTION:  1. 


Set  hairline  of  cursor  on  Q  value  of  1000  on  D  scale. 


2.  Move  center  slide  until  assumed  stability  parameter,  n, 
of  0.25  on  cursor  intercepts  assumed  wind  speed,  IT,  of 
2  m/sec  on  C  scale, 

3.  Move  cursor  until  n  value  of  0.25  on  cursor  intercepts 

mg  min 

desired  dosage  of  100  — ^  on  A  scale. 

°  cum 

4.  Read  area  on  B  scale  indicated  by  hs.irline  of  cursor 
(A  =  610  sq  m), 

IV.  CONCLUSIONS. 

Slide  rules  have  been  contrived  for  evaluating  O.  G.  Sutton's 
mathematical  models  for  estimating  the  downwind  travel  and  area  coverage 
from  an  instantaneous  point  source  of  vapor  or  very  fine  aerosol  disseminated 
and  sampled  at  ground  level. 

Using  these  slide  rules,  the  nonmathematician  can  employ  the 
mathematical  theory  for  making  preliminary  weapons -effects  predictions 
necessary  in  handling  everyday  problems  encountered  in  research  and 
development  of  chemical-weapons  systems. 
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